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Distinguishing black-hole spin-orbit resonances by their gravitational wave signatures. 

II: Full parameter estimation 
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Gravitational waves from coalescing binary black holes encode the evolution of their spins prior 
to merger. In the post-Newtonian regime and on the precession timescale, this evolution has one of 
three morphologies, with the spins either librating around one of two fixed points (“resonances”) or 
circulating freely. In this work we perform full parameter estimation on resonant binaries with fixed 
masses and spin magnitudes, changing three parameters: a conserved “projected effective spin” £ 
and resonant family A4> = 0, tt (which uniquely label the source); the inclination Ojn of the binary’s 
total angular momentum with respect to the line of sight (which determines the strength of preces- 
sional effects in the waveform); and the signal amplitude. We demonstrate that resonances can be 
distinguished for a wide range of binaries, except for highly symmetric configurations where preces- 
sional effects are suppressed. Motivated by new insight into double-spin evolution, we introduce new 
variables to characterize precessing black hole binaries which naturally reflects the timescale sepa¬ 
ration of the system and therefore better encode the dynamical information carried by gravitational 
waves. 

PACS numbers: 04.25.dg, 04.70.Bw, 04.30.-w 

I. INTRODUCTION 

Gravitational waves (GWs) from binary black holes 
(BBHs) are expected to be an important source [IH4] 
for future networks of GW detectors such as Advanced 
LIGO/Virgo (5][6], LIGO-India [7], KAGRA 00, and 
the Einstein Telescope m- Once formed, BBHs emit 
GWs that extract energy and angular momentum from 
the orbit, decreasing the binary separation and increas¬ 
ing the orbital frequency (and thus the GW frequency). 

Most binaries are expected to circularize by the time 
they enter the sensitivity band of ground-based detec¬ 
tors mm (see e.g. HM7] and references therein for 
recent work on eccentric binary waveforms, rates and 
detection strategies). Circular BBH inspirals are char¬ 
acterized by the masses (mi, m 2 ) and spins (Si, S 2 ) of 
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each black hole. The spin and orbital angular momenta 
precess on timescales shorter than the radiation-reaction 
timescale [T8H2Q] . Their general evolution can be un¬ 
derstood relative to the two one-parameter families of 
fixed points (“post-Newtonian resonances”) of the orbit- 
averaged spin precession equations [21]. At resonance, 
the spins and orbital angular momentum L remain copla- 
nar. If A4> denotes the angle between the projection 
of the spins in the orbital plane orthogonal to L, the 
spins can either lie on the same side of L (the A<f> = 0 
resonance) or on opposite sides of L (the A4> = 7 r reso¬ 
nance). In general, the conservative evolution of the spins 
of any BBH can be understood as either a trapped phase- 
space orbit around one of the two fixed points (A<F = 0 
and A^> = 7 r) or as an orbit where A4> circulates freely. 
Therefore the phase space can be split into three classes, 
or “morphologies” [221123] , 

The morphology of BBHs in the Advanced 
LIGO/Virgo band is determined by the spin con¬ 
figuration when the compact binary is first formed, long 
prior to its detection via GWs [24] . This relationship 
may enable GW measurements to constrain the efficiency 
of tidal interactions in binary star evolution and may 
indicate whether a mechanism such as mass transfer, 
stellar winds, or supernovae can induce mass-ratio 
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reversal (so that the heavier black hole is produced by 
the initially lighter stellar progenitor). More broadly, 
spin measurements of BBHs have long been expected 
to be a critical element in distinguishing the underlying 
astrophysical model for compact binary formation, 
complementing other measurements: see e.g. [251128] . 

Gravitational radiation encodes the dynamics and 
properties of the inspiralling binary. Precessing bina¬ 
ries produce a rich, highly modulated signal mmana- 
mi, enabling constraints on compact binary parameters, 
such as masses and the misaligned spins [35H4T] . This in¬ 
formation can be recovered from detector measurements 
by systematically comparing all possible candidate sig¬ 
nals with the data and constructing a Bayesian posterior 
probability distribution for all binary parameters |4lT[48] . 
Out of all precessing BBH configurations possible a pri¬ 
ori, the two librating morphologies are relatively common 
for comparable-mass binaries in the close orbits to which 
LIGO is sensitive, but relatively rare for binaries drawn 
from a uniform mass distribution 123] . 


That said, because the relative spin dynamics do not 
directly source the dominant GW signal, it is important 
to investigate the extent to which the relative orientation 
of precessing spins can be constrained in general, and the 
degree to which measurements distinguish morphologies 
in particular. Vitale et al. m recently studied this 
problem for a small sample of precessing binaries, claim¬ 
ing that the relative angle A4> between the projection of 
the two spins into the orbital plane cannot be measured. 


Previously, Gerosa et al. [49] examined the GW signal 
from the two resonant families and evaluated the over¬ 
lap between the waveforms with all parameters fixed, to 
demonstrate that the two families produce qualitatively 
and quantitatively distinct GW signatures. Gupta and 
Gopakumar [501 used overlaps to imply resonant and non¬ 
resonant binaries were distinguishable. In this work we 
apply the state-of-the-art lalinference parameter es¬ 
timation code [48] to the expected detector response due 
to resonant post-Newtonian binaries in the nearby Uni¬ 
verse. Using the posterior distribution of compact binary 
parameters, we assess for the first time how confidently 
these measurements distinguish between the three possi¬ 
ble morphologies. As we discuss in detail, morphological 
classification is related to (but not strictly dependent on) 
the accuracy with which we can measure the angle A<f> 
between the projection of the two spins on the orbital 
plane. On physical grounds, morphological classification 
is far more robust: while A4> changes on the precession 
timescale, a binary’s morphology only changes on the in¬ 
spiral timescale. 

The paper is organized as follows. In Section [IT] we 
briefly review GW parameter estimation, the specific grid 
of binaries examined for this work and the parameters 
needed to interpret the posterior distribution in the con¬ 
text of post-Newtonian resonances. In Section [HI] we 


show that the morphology of exactly resonant binaries 
can be reliably determined at astrophysically plausible 
signal amplitudes, unless the binary is face-on. Moti¬ 


vated by a recent analytic solution of the spin preces¬ 
sion equations [22] [23] , we introduce two coordinate sys¬ 
tems to more naturally characterize how well properties 
of double-spin binaries have been constrained. We con¬ 
clude in Section [IV] describing how our results impact 
the broader program of astrophysical inference using GW 
measurements. In an attempt to keep this paper self- 
contained, we report various technical material in the ap¬ 
pendices. Appendix [A] describes the relationship between 
several conventions used to describe precessing spins that 
have appeared in the literature. Appendix [B] reviews 
the characteristic precession timescales introduced in our 
previous work on analytic solution of the two-spin pre¬ 
cession equations [22l 123], Finally, in Appendix |C| we 
provide some technical details on our parameter estima¬ 
tion procedure. 


II. METHODS 

A. Candidate sources and signals 

The GW signal from the binary depends on 8 intrinsic 
(physical) parameters - the component masses mi and 
m 2 and spin vectors Si, S 2 (3 components each) - and 7 
extrinsic parameters - the event time £ re f; luminosity dis¬ 
tance Dl ; sky location or propagation direction relative 
to the Earth’s equatorial coordinate system, expressed 
via two angles <a, S (right ascension and declination); and 
three Euler angles describing the orientation of the bi¬ 
nary relative to the Earth. The rotation relating the 
Earth frame and the binary’s frame is commonly char¬ 
acterized using the line of sight N from the observer to 
the binary and the time-dependent orbital angular mo¬ 
mentum direction L, evaluated at some fiducial orbital 
frequency. For example, the orbital orientation is of¬ 
ten characterized by the inclination angle defined by 
cos l = — L • N (in this convention, a binary whose angu¬ 
lar momentum points towards the observer corresponds 
to l = 0.) Among different possible parametrizations for 
these intrinsic and extrinsic parameters, we choose to use 
the following [55] : 


6 = {m!, m 2 , xi, y 2 , 0 JN , 0 LSl , 0 LS2 , A4>, </> JL } , (1) 

as well as coalescence event’s 4 spacetime coordinates, 
polarization, and orbital phase, which will be marginal¬ 
ized over and not discussed henceforth. In this list, the 
symbol m^ denotes the component masses, Xi — 
the adimensional spin magnitudes, OlSi an gl e be¬ 
tween the orbital angular momentum L and the spin vec¬ 
tor Si, Ojn the angle between the total angular momen¬ 
tum J = L + S 1 + S 2 and the line of sight to the observer 
(—N), A4> the azimuthal angle between Si and S 2 , and 
(j)jL the azimuthal angle between J and L, relative to 
a reference phase defined via the projection of —N onto 
the orbital plane; see Figure [T] and 39.. All parameters 
are specified when the GW frequency (twice the orbital 
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frequency) is / re f = 100Hz. Except for mi, m 2 , xi? X2> all 
of our intrinsic parameters evolve during the inspiral. 

The orbital dynamics and GW signal are constructed 
via the lalsimulation SpinTaylorT4 code [ 51 ] . Based 
on previous implementations [521 [53], this time-domain 
code solves the orbital dynamics of an adiabatic, quasi¬ 
circular inspiraling binary using the “TaylorT4” method 
[54] for the phase evolution and (orbit-averaged) pre¬ 
cession equations for the angular momenta [19]. In all 
cases we used a GW phase up to 3.5PN in nonspin¬ 
ning terms; up to 2.5PN in spin-orbit terms; up to 2PN 
spin-spin terms; and orbit-averaged spin precession equa¬ 
tions up to 2PN, including the leading-order spin-orbit 
and spin-spin interactions. As described below, we per¬ 
formed calculations both with and without black hole- 
consistent coefficients in composition-dependent terms 
like the quadrupole-monopole (QM) coupling. We evalu¬ 
ate the GW amplitude h{t) using only the leading-order 
(Newtonian) quadrupole [55] . 

To compare waveforms from binaries belonging to dif¬ 
ferent morphologies in controlled conditions we inves¬ 
tigate only a specific combination of the binaries’ in¬ 
trinsic parameters: total mass M = 13.5 Mq, mass ra¬ 
tio q = m 2 /mi = 0.8 (which implies m 2 = 6.0M 0 
and mi = 7.5M 0 ), and adimensional spin magnitudes 
Xi = X 2 = 1 (maximally spinning BHs). Our choice, 
motivated by population-synthesis predictions [HU and 
computational restriction^] also facilitates comparisons 
with our previous works [24[ [49 . 

Following [49] , we parametrize binaries from each mor¬ 
phology using the dimensionless projected effective spin 
£ 


= S 0 • L = Xi cos #LSi + Q X 2 cos 6 LS2 ( , 

~ M 1 2 f=f rei 1 + q 

where So is the effective spin 

So = (1 + q) Si + ^1 + S 2 . (3) 

The projected effective spin £, which was first introduced 
in [561157] and used to solve the spin precession equations 
in E 2 U 22 ], is a constant of motion at 2PN order in spin 
precession and 2.5PN order in radiation reaction. 

Our candidate sources will be exactly resonant, which 
is a conservative choice, since in this case the waveforms 
are weakly modulated and our ability to constrain param¬ 
eters correlates with the modulation in the waveforms. 
For resonant binaries, for each £ and each family there are 
unique values of 6lSd @ls 2 i AT, meaning that (£,AT) 
uniquely label the source. For simplicity, we fix c/)jl = 0 
at / r ef = 100Hz, so N,L,J are coplanar. Finally, be¬ 
cause the amount of spin-induced precession seen by the 


1 Data presented in this paper involves full MCMC simulations 

which required several million CPU hours over the course of sev¬ 
eral months to compute. 




FIG. 1. Definitions of the angles used throughout this paper. 
The directions of the spins Si and S 2 are specified using polar 
angles OlSi and azimuthal angles T* (i m 1,2), measured 
relative to x ; = x x L. For resonant binaries, the orbital 
angular momentum and spins remain coplanar, implying that 
the angle AT = <f >2 — Ti remains constant: either A<f> = 0 or 
AT = 7r. 
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Parameter 

Value (s) 

mi 

7.5 M© 

m 2 

6 Mq 

Xi 

1.0 

X 2 

1.0 

L/M 2 

0.896 

A<f> 

0, 7T 

£ 

0,JN 

-1, -0.9, ..., 1 

A 7T f 37T 7T 

U ’ 8 ’ 4 ’ 8 ’ 2 


TABLE I. Overview of simulated binaries. Input binary 
parameters at / = 100Hz (equivalent to v = y/nfM — 0.208, 
where A4 = (m\rri 2 ) 3 ^/(rrti + m 2 ) 1//5 is the chirp mass) for 
the 2 x 21 x 5 distinct parameter estimation calculations de¬ 
scribed in the text. 


observer in the GW signal depends on Ojn [49] . we ex¬ 
plore a range of Ojn between 0 and tt/2 without loss 
of generality: the posterior distribution for sources with 
0jn > 7t/2 can be related to a posterior with Ojn < tt/2 
by suitable reflection symmetry. Hence, having chosen 
the resonant family (A4> = 0 or A4> = 7r); the inclina¬ 
tion angle Ojn (0, i r/8, 7t/4, 37t/8, 7t/2); and £ in the 
[—1,1] interval in steps of 0.1, we obtain 42 candidate 
sources and 210 candidate signals. 

For simplicity we adopt the sky location (<a, S) = 
(0°,0°) and fiducial GPS (Global Positioning System) 
time equal to zero. To better isolate the intrinsic differ¬ 
ences between these signals, we require all sources to pro¬ 
duce the same coherent amplitude in a network consisting 
of the three (two Advanced LIGO plus Virgo) interfer¬ 
ometers, using the design sensitivity curves [58] . Specifi¬ 
cally, we select the distance for each source such that the 
network signal-to-noise ratio (SNR) is 20, a high but not 
unreasonable value for first detections HHSiEo]. 

This procedure defines two families of candidate 
sources, one for A4> = 0 and one for A4> = 7r, each 
of which is characterized by two parameters (£, defined 
in Eq. [ 2 j and 6jn), and produces a well-defined strain 
response in each idealized instrument. For all masses 
and spins, the family with A4> = i r has a special point 
£ = £t = (xi ~ <?X2)/(l + (z) such that both spins are par¬ 
allel to the orbital angular momentum and antialigned 
with respect to one another, i.e. ( 0ls 1 ^ls 2 ) = (TP 0 ) 
EQ- For our choice of binary parameters, £^ —0.11 

[see Eq. (J2| and Figure [2]. 

Rather than explore an ensemble of noise realizations, 
following previous studies dunums we adopt a unique 
preferred noise realization (exactly zero noise) for all sim¬ 
ulations and all instruments; see Appendix [C] for details. 
As discussed below, we have repeated parts of the analy¬ 
sis using randomly selected detector noise to demonstrate 
that our results are robust with respect to this choice. 
Simulation properties are summarized briefly in Table [T| 
in the Appendix, Table [TT] provides an explicit list, includ¬ 
ing derived parameters (such as J and the termination 
frequency). 

As reviewed at length in [49] , each one-parameter fam¬ 
ily of resonances has distinctly different angular momenta 


orientations as a function of £: in other words, varying 
£ significantly changes the dynamics of the source bi¬ 
nary. For example, as shown in the left panel of Figure 
[2] our one-parameter family of dynamically unique bina¬ 
ries with A<f> = 0 has both spins comparably (costs'! — 
cos 0ls 2 ) and often significantly misaligned with L. As 
a result, the A4> = 0 binaries have significant misalign¬ 
ment between L and J, as seen in the right panel of 
Figure |2j By contrast, our family of A4> = n sources all 
have either one or the other spin nearly aligned with the 
orbital angular momentum: cos 0ls 2 — 1 for £ > £t and 
cos 6ls\ — — 1 for £ < As shown in the right panel of 
Fig-i our A4> = 7r binaries have minimal misalignment 
between L and J, particularly when £ < £t- We will see 
below that some measurements are strongly correlated 
with the amplitude of precession-induced modulations, 
and therefore correlate strongly with the misalignment 
angle #jl(£) = arccos J • L. Note also that, for our sim¬ 
ulations, the values £ = ±1 correspond to spin-aligned, 
non-precessing binaries. 


The quadrupole-monopole term 


Consistent with all previous parameter estimation 
studies and the implementation of lalsimulation at 
the time this project began, our spin precession equa¬ 
tions explicitly omitted the QM coupling described in 
m- After realizing the omission, we have added the 
missing terms to the lalsimulation code, and we have 
repeated parts of our computationally expensive analysis 
including the QM interaction. The conclusions of this 
exercise are encouraging, and our main results are very 


robust under perturbations: as discussed in Sec. IIIF be¬ 
low, the posterior distributions obtained with and with¬ 
out the QM interaction are basically indistinguishable, 
at least for near-resonant binaries. 

The omission of the QM term has several important 
consequences from a theoretical point of view: (i) The 
projected effective spin £ is only approximately conserved 
in our simulated binaries. If £ is not conserved, spin pre¬ 
cession in general will not be quasi-periodic. Fortunately, 
as argued in Sec. EH below, £ is nearly conserved even 
in the absence of the QM interaction for near-resonant 
binaries, so this does not have a dramatic impact on the 
present analysis, but the omission of the QM interac¬ 
tion could have larger effects for binaries that are far 
from resonance. To eliminate any ambiguity, £ (like all 
spin-dependent intrinsic parameters) must be evaluated 
at some reference frequency, here chosen to be 100Hz. (ii) 
The morphological classification described in j22[ [23] is 
no longer exact: for example, binaries that should be ex¬ 
actly resonant can, in practice, perform small librations 
around resonant configurations, and some binaries may 
not belong to any of the three morphologies discussed 
in (22JI23]. (iii) Precessional dynamics is obviously mod¬ 
ified at 2PN, and so is the evolution of the binary on the 
radiation-reaction timescale. 
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FIG. 2. Review of resonant dynamics and our simulations. Left panel: Illustration of the specific two one-parameter 
families of resonant binaries used in this work. The two curved paths show cos OlSx , cos 0ls 2 for A<I> = 0 (blue and dotted) 
and = n (red and solid), superimposed on contours of constant £ (dashed). £ = £ T — —0.11 is shown with a thick dashed 
line. Stars mark the location of each of our simulations (blue for A<f> = 0, red for A<I> = 7r). Right panel : Plot of J • L versus £ 
for each resonant family at the reference frequency / re f = 100 Hz. As discussed in the text, our source binaries in the A<1> = n 
resonance have spins and L much more closely aligned with each other than sources in the A<1> = 0 resonance. All source 
binaries have L > Si + S 2 , and hence dynamics characterized by the two-spin solution described in [22l [23] , 


B. Gravitational wave parameter estimation for 
double-spin binaries 

Parameter estimation can be performed using Bayes’ 
theorem: a prior probability distribution p{0\H) (where 
H is the hypothesis and 6 represents collectively the pa¬ 
rameters of the system) is updated upon receiving data 
d from the experiment to give a posterior distribution 
p(6\d,H) 


p(9\d,H) = 


gMjP) 

p(d\H) 


(4) 


where p(d\Q, 'LL) is called the “likelihood”. Starting from 
an appropriate prior distribution, samples of the poste¬ 
rior distribution are randomly selected through a stochas¬ 
tic sampler using information from the data. We per¬ 
form parameter estimation on simulated GW signals 
from resonant binaries using the lalinference suite 
[48], in the lalinference_mcmc implementation. This 
Markov Chain Monte Carlo (MCMC) algorithm gener¬ 
ates a series of independent samples from the posterior 
distribution, given prior probabilities for all physical pa¬ 
rameters. Consistent with previous work [48] . we a priori 
assume component masses are uniformly distributed be¬ 
tween 1 Mq and 30M 0 , with total mass less than 35M©; 


component spins randomly oriented, with uniformly dis¬ 
tributed amplitude; and the merger event occurs within 
±0.1s of our proposed event time t re f. To account for the 
high mass of our candidate sources and the high sensi¬ 
tivity of the advanced instruments used in this study, we 
allow the physical position to be uniformly distributed 
on a sphere of Euclidean radius 6Gpc, implying unifor¬ 
mity over the sky and in distance - i.e., p(d) oc d 2 - out 
to 6Gpc. The prior was not modified to model selec¬ 
tion bias. This prior includes extremely high mass ratios 
and a wide range of spins, favoring binaries dominated 
by a single spin rather than systems showing strong two- 
spin effects; our conclusions drawn using this unfavorable 
and astrophysically unmotivated prior should therefore 
be taken as conservative underestimates. For simplic¬ 
ity and consistent with previous work, cosmological ef¬ 
fects are not included in our calculations^] Finally, we 
stop our calculations when the final Markov Chain, after 


2 The most distant sources in our sample, with #jiv — 0 and 
£ ~ ±1, were placed at cT ~ 815Mpc, or 2 : ~ 0.2 with current 
cosmological parameters; typical sources have cIl — 500Mpc, or 
2 : ~ 0.12. Because the GW strain depends only on the redshifted 
masses, the posterior including cosmological effects can be de¬ 
rived from a cosmology-free posterior by a coordinate transfer- 
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FIG. 3. The J£-plane. For fixed source parameters rai, 7712 , Xi> X 2 and / re f, the solid black lines indicate the extent of the 
allowed values of J and £. Any point on the lower (upper) edge of this region corresponds to a A<f> = 0 (A<f> = ir) post- 
Newtonian resonance. The specific parameter configurations used in this paper are illustrated by red stars (on the upper 
boundary) and blue stars (on the lower boundary), as in Figure [ 2 ] Left: For these same fixed source parameters, three shaded 
regions (red on top; blue on the bottom; and green in between) show the regions for, respectively, librating about A<f> = 7 r, 
about A<f> = 0, and circulating precession morphologies (described in the text). Dashed lines separating these regions (marked 


by arrows) show morphology boundaries JcO,ctt(£)- Right: Prior probability p(J, £), derived from the priors described in IIB 


by marginalizing out mi, m 2 , xi, X 2 - For scale, the black curve also shows the same boundary of allowed (J, £) values shown in 
the left panel. The structure in this distribution reflects the broad range of mass ratios and spin orientations included in the 
prior; see the text for details. 


downsampling to remove correlations, has roughly 2,000 
uncorrelated samples from the posterior [48] . 

The parameters of the posterior samples can be repre¬ 
sented in any coordinate system, allowing us to efficiently 
compute the posterior p(Q) using any parameters 0 , inde¬ 
pendent of the coordinates used internally by the MCMC 
simulation. The relative orientation of the two spins with 
respect to the orbital angular momentum, OlSx and @ls 2 5 
can be used to characterize the system. Motivated by 
the analysis of [22] [23] and as described at greater length 
in Appendix [A] we replace these angular variables with 
(J, £), where £ is given by Eq. © and 

J = [L 2 + S\ + Si + LS\ cos e LSl + LS 2 cos 0 L s 2 
+ 2 S 1 S 2 (sin 0 LSl sin 0 L s 2 cos A< 1 > + cos 0 LSl cos Q L s 2 ) ] 5 

(5) 

is the magnitude of the total angular momentum, J = 
|J|. For each resonant family (A<£> = 0,7r), each pair 


mation, using a modified (Euclidean) distance prior rather than 
one consistent with cosmology. At the distances of interest, a 
Euclidean and cosmological distance prior nearly agree. 


(J, £) uniquely specifies the binary configuration (see 
Fig.§. The magnitude L = |L| of the orbital angular 
momentum is calculated at leading (Newtonian) order: 
L = r]M 2 /v = 77M 2 /(7r/ re fM) 1 / 3 . Unlike the system- 
frame parameters described above, these parameters nat¬ 
urally reflect the separation of timescales in the two- 
spin problem, with £ conserved up to 2PN order on all 
timescales by the orbit-averaged spin-precession equa¬ 
tions; J changing on the radiation-reaction timescale; 
and A<£> changing on the precession timescale [22] [23] . 
For each fixed mi, 777 , 2 , Xi> X 2 and /ref, a range of (J, £) 
is allowed; see Gerosa et al. [23] for details. For single¬ 
spin binaries, the relationship between J and £ at fixed 
L is one-to-one; for double-spin binaries, a range of J 
are allowed at each fixed £. To guide the eye in the 
plots that follow, we evaluate and show this “allowed” 
region in the J£-plane for the chosen source parameters 
mi,m 2 ,Xi,X 2 in Figure^ We will use J r0 (0,-WO to 
denote the maximum and minimum values of J for a 
given £. Points on these two curves, including the injec¬ 
tions described above, are all resonant binaries: binaries 
with J = J r7T (<£) belong to the A<£> = i r resonance, while 
those with J = J r o(£) belong to the A4> = 0 resonance. 

Additionally, as described in Refs. [22] and [23] . these 
parameters facilitate morphological classification , sub- 
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dividing the parameter space H into three disjoint re¬ 
gions HcnHc,^, set by the (nondissipative) dynamics 
at / re f = 100Hz. Specifically, Ho and H^ are the systems 
whose spins are librating about the A4> = 0 and A4> = 7 r 
resonance, respectively; He are the remaining, circulat¬ 
ing binaries. Geometrically, these three regions are sepa¬ 
rated by configurations where either Si or S 2 are parallel 
to L at some point during one precession cycle. From a 
computational point of view, motivated by [22j [23] , at 
each mi,m 2 ,Xi,X 2 we determine the morphology from 
the values of J and £ at 100 Hz. We define J c o(£), T C7r (£) 
as the two values of J which allow either Si or S 2 to 
be parallel to L at some point during their precession 
cycle, ordered as J r7r (£) < J C7r (£) < J c o(£) < T r0 (£). Be " 
cause these surfaces separate the different morphologies, 
we identify the appropriate morphology simply by com¬ 
paring J to JcO,cvr(0 : binaries with J < J C7r (£) belong to 
H^; binaries with J > J c o(£) belong to Ho; and binaries 
with J C7r < J < J c 0 belong to He- These boundaries are 
illustrated in Figure [3j 

We apply morphological decomposition to the list of 
posterior samples, classifying the fraction p(A) of sam¬ 
ples that are consistent with each morphology A = Ho, 
Htt, and He, respectively. As described in previous work 
(see in particular [2 23]), while librating binaries occur 
frequently near merger for comparable-mass binaries, a 
relatively small amount of the possible parameter space 
describing all compact binaries corresponds to librating 
orbits trapped near a resonance (H^o), as opposed to 
circulating orbits (He)- 

Usingthe prior parameter distribution described below 
Eq. fflrlwe find prior probabilities p pr ior(Ho, 7 r,c) 

Pprior('Ho) = 0 . 118 ; 

Pprior(74) = 0.049; (6) 

Pprior('Hc) = 0.832 . 

The prior also implies a marginal distribution p(J, £), 
shown in the right panel of Figure [3] The marginal distri¬ 
bution has several salient features, most notably a peak 
near (J, £) ~ (0.6,0) arising for example from randomly- 
oriented spins in unequal-mass binary black holes. Ex¬ 
tending away from that peak are three wedges, corre¬ 
sponding to binaries (a) with spins increasingly aligned 
with L, extending upward toward J ^ 1 and £ ~ 1; (b) 
with spins increasingly anti-aligned with L, extending 
downward towards £ — 1; and (c) a region extending 


3 To ensure consistency, we draw samples from the prior by run¬ 
ning the parameter estimation code with input h{t) = 0. While 
we are forced to adopt the ad-hoc choices of Eq. © for prior 
probabilities to be consistent with the mass and mass ratio dis¬ 
tribution described previously and used by convention in lalin- 
ference, we emphasize that those probability distributions were 
chosen arbitrarily to avoid the impression of bias in mass. Nature 
likely favors other mass distributions mu, which may further 
enhance the prior probability that non-circulating morphologies 
occur. 


to the left along £ ~ 0 with small J, corresponding for 
example to small mass ratios q and small primary spins 
Xi- 

The posterior distribution can also be parametrized 
using the two frequencies implied by the spin precession 
equations, evaluated (to avoid ambiguity) using a ref¬ 
erence orbital frequency: here / re f = 2/ or b = 100Hz. 
These two timescales appear directly in the spin and or¬ 
bital dynamics, and hence are reflected in the waveform. 
Rather than adopt a precession-phase dependent defini¬ 
tion of these two timescales via, e.g., the vector that 
appear in the spin precession equations dSi/dt = x S$, 
following Kesden et al. [22] . Gerosa et al. [23! we use 
an analytic, explicitly doubly periodic solution for the 
nondissipative precession dynamics to identify a preces¬ 
sion timescale r for the relative spin motion [defined be¬ 
low Eq. (8) in [251] and a precession frequency (Q z ) for 
the motion of L around J. Relative to their notation, 
(Q z ) = a/ r, where a is the precession angle of L around 
J in one spin precession period r; see Appendix [B] for 
details. 


III. RESULTS 

A. Robust morphological classification 

Let us first consider sources with A4> = 0. Except 
for binaries that are either face-on (6jn = 0) or have 
both spins nearly aligned with L (| cos OlSi \ — 1 and 
| cos 0ls 2 I — 1 ), oar calculations show that the posterior 
probability p(Ho) is a good indicator of which resonant 
family the source came from, for almost all source orien¬ 
tations and £, and for astrophysically plausible source 
amplitudes. For example, the left panel of Figure [4] 
shows the posterior probability p(Ho|A4> = 0 ,£,6 >jat): 
the fraction of posterior probability implying the source 
librates about the A<F = 0 resonance, as a function of 
the true source parameters (£,^jat) at a fixed source 
amplitude SNR = 20. This probability is large (and 
nearly unity) for almost all (£, 0 jjv)- For binaries at or 
near the A<F = 0 resonance, we can confidently classify 
them in their corresponding resonant regime. Equally 
striking, the center-left panel shows that the probability 
of misidentification in the opposite resonance is nearly 
zero: p(H 7T \A <F = 0, £, 0jjv) — 0- These high prob¬ 
abilities should be sharply contrasted with the small 
prior probability associated with the A<F = 0 resonance 
[Pprior(Ho) = 0.118], as the majority of phase space is 
associated with circulating binaries. 

The right panels of Figure [4] show our corresponding 
but qualitatively different results using source binaries 
from the A4> = 7 r resonant family. Except for specific 
configurations, our calculations show that the posterior 
probability favors the circulating morphology and dis¬ 
favors the “incorrect” classification. This follows from 
geometrical considerations: A4> = 7 r sources are more 
aligned and precess less, making them more difficult to 



A<3> = 7 r sources 



£ 


£ 


FIG. 4. Constraining morphology versus inclination: the posterior. Classification probability versus £ for different 
values of Ojn at fixed SNR = 20. Panels in this plot are organized as a matrix: sources (A<f> = 0 or A<f> = ty) are columns; 
the fractions of the posterior in each morphology are rows; and each panel shows the corresponding probability 

evaluated for sources at different values of £. For example, the bottom-left panel shows the probability p(7ic |A<f> = 0) versus £ 
and (in different line styles) versus Ojn- Nearly aligned binaries (£ ~ 1,-1 or A<f> = ty and £ = £t — —0.11) and nearly face-on 
binaries (Ojn — 0) cannot be reliably classified (top left and middle right). A source with A<1> = 0 can be reliably classified 
in the Tio morphology (top left) unless face-on (for our simulations, £ « ±1); however, a source with A4> — ty can be reliably 
classified in the morphology only if Ojn — ty/ 2 and £ ^ — 1,1, £t (green curve in middle right). In other words, cases where 
the probabilities vary significantly with £ reflect how the misalignment angles between the various angular momenta change 
versus £: compare to Figure [5] Due to the relatively small number (N ~ 2000) of posterior samples used, these probabilities 
have a few percent (binomially distributed) sampling uncertainty. 
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FIG. 5. Constraining A<I> at 100 Hz. Probability distribu¬ 
tion dP/dA$> of the phase angle A<f> between the spins can 
be reliably measured for A<f> ~ 0. Plot shows the posterior 
distribution for A<I> between the spins using spin vectors at 
the reference frequency / re f = 100 Hz, given data containing 
sources with A<f> = 0 (left column) or A<1> = tt (right column), 
with either Ojn — 7t/2 (top) or 6jn — 7t/4 (bottom). In all 
panels, colors indicate different values of £, as shown in the 
legend. This figure shows that the posterior distribution in 
the intrinsic parameter A<f> depends strongly on the extrin¬ 
sic parameter Ojn and weakly on the intrinsic parameter £, 
except for nearly aligned systems £ ~ ±1. 


distinguish from non-precessing sources. The probability 
of “correct” classification p(J~L n \ A= tt, £, Ojn ) shown in 
the middle-right panel decreases significantly at £ = £t — 
—0.11. This value of £ corresponds to a special configu¬ 
ration featured by A<4> = i r sources, with cos^Si = —1 
and cos 0ls 2 — lj see, e.g., our Figure [2] or Figure 2 of 
[49] . For our injection masses and spin magnitudes, this 
“aligned” configuration occurs at £ = £t — —0.11, as ex¬ 
plained earlier^] Resonant binaries at and near £^ have 
spins that barely precess, being almost parallel to L. In 
other words, just like the two limits £ -A 1 and £—>►—!, 


4 The spin-aligned configuration at £ = on the A<f> = 7r family 
arises because the down/up configuration is a stable fixed point 
of the spin precession equations; by contrast, for our masses, 
spin magnitudes, and frequencies, the up/down configuration is 
unstable and not a member of the A<f> = 0 family j61| . When 
the up/down configuration is stable, the point with cos$i = +1 
and cos 02 = — 1 will belong to the A4> = 0 resonant family and 
we anticipate similarly diminished precession amplitudes in its 
vicinity leading to a similarly difficult challenge in identifying 
the morphology. 


when both spins are coaligned with each other and L, 
resonant sources with A<4> = n and £ £ T have almost 
indistinguishable dynamics from a non-precessing binary. 
As seen in Figure [4] for sources with £ ~ £t and A<4> = tt 
the morphology is only weakly constrained by observa¬ 
tion: p(H 7r|A<4> = 7T, —0.11 , Ojn) is small, comparable 
to the prior (see Sec. Ill E| and Fig. [I2| below). Second, 
the probability for correct classification of a A<£> = i r 
source is often significantly smaller than the correspond¬ 
ing A<£> = 0 result. As emphasized in [49], the A<£> = i r 
and A<£> = 0 resonant families have qualitatively different 
dynamics and GW signals: A<4> = 7r resonances undergo 
substantially less precession of L, and hence have less 
modulated GWs. As a result, at fixed signal amplitude, 
GW measurements of sources from the A<4> = i r reso¬ 
nance should be less effective at measuring double-spin 
physics like the morphology, implying the lower probabil¬ 
ities p(H 7 r |A<F = 7T, £, Ojn) < p(Uo |A$ = 0,£,6> JA r) seen 
in the corresponding panels of Figure [4] 

To qualitatively corroborate this analysis, we also ex¬ 
amined the posterior for A<4>; an example appears in Fig¬ 
ure [5] In most cases, A<4> could be weakly constrained to 
be 0 or ~ 7T. By definition, the region 1~Lq corresponds 
to libration about the A<4> = 0 resonance. Because the 
phase-space trajectories in this region explicitly cannot 
extend to A<£> = tt [22] [23], posterior samples in this re¬ 
gion explicitly disfavor A<£> tt; in fact, depending on 
£, the allowed region of A<£> can be narrowj^] To summa¬ 
rize, our ability to confidently classify source morphol¬ 
ogy to Ho, 7 r implies A<£> can be constrained away from 
7T, 0, respectively. That said, morphological classification 
is more robust and qualitatively different than measure¬ 
ments of A<4>; for example, while A<4> changes on the pre¬ 
cession timescale, a binary’s morphology only changes on 
the inspiral timescale. 

All of our results depend on the signal amplitude, here 
fixed at SNR = 20. Repeating parts of our study with 
much stronger signals (SNR ~ 50), we have found dra¬ 
matically improved morphological classification, even for 
sources with A<4> = 7r. Conversely, our ability to reliably 
classify morphologies gradually degrades at lower source 
amplitude, as shown in Figure [6] 


B. Using conserved quantities as coordinates to 
describe double-spin measurements 

On physical grounds, we expect GW measurements to 
best constrain quantities that are constant on the longest 


5 At sufficiently high signal amplitudes, phase angles like A<f> can 
be measured directly, independent of how rapidly they evolve 
with frequency. For example, the phase angle </)jl can be es¬ 
timated to within an angle of order 27r/SNR, with a coefficient 
that depends on the line of sight (Ojn) and the amplitude of 
precession of L around J (e.g., the characteristic angle between 
L and J): see, e.g., m- 
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£ £ 

FIG. 6. Constraining morphology versus SNR. Plot of the classification probability versus £ for different choices of SNR 
and Ojn — 7t/4. The panels in this figure are organized as in Figure [4] As described in Appendix [C] the SNR = 20 results (red 
curve) were produced using a lower starting frequency than results for SNR< 20; this discrepancy is responsible for the slight 
difference in trend between the right panels’ results for SNR=20 versus results for SNR< 20. 
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FIG. 7. Constraining (J, £). The parameters J and £ can be independently constrained to a narrow volume of parameter 
space. Posterior distributions from distinct simulations in the J£-plane for 6jn — tt/ 4 and SNR = 20, for source binaries with 
A<f> = 0 (left) and = i r (right). The actual value from each simulation is marked with a star; the maximum-likelihood 
estimate is marked by a black cross. To guide the eye, the allowed region for J and £ for the true source parameters is 
superimposed (thick black line) as in Figure [3] As the allowed region depends on masses and spins, not all points in the 
posterior lie in this region; see i23 for details. Solid and dashed black lines show the 95% and 67% confidence intervals for each 
set of posterior samples from distinct simulations. Top panels: Samples are colored according to morphology (blue: librating 
around A<f> = 0, red: librating around A<£> = 7r, green: circulating). Bottom panels: Samples are colored according to the 
posterior mass ratio q. In our simulations, binaries without sufficiently strong precession-induced modulation to enable tight 
constraints on q are equally poorly constrained in morphology: our ability to classify morphology correlates with our ability to 
constrain q. 
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timescales. In a two-spin system, while the specific com¬ 
ponent spins and orbital angular momenta depend ef¬ 
fectively on a precession phase choice, recent analytic 
work [22] [23] has identified natural quantities, conserved 
on the precession timescale, to characterize the orbit: 

(j, o- 

Figure [5] shows posterior probability distributions for 
A<h, for several fiducial source events; Figure [7] shows 
corresponding posterior distributions for (J, £). To guide 
the eye, in Figure [7] a thick black line shows the allowed 
region for J, £ assuming all mass and spin parameters are 
equal to the source parameters, as in Figure [3] First and 
foremost, as noted above, the measurements shown in 
Figure [5] immediately reveal that A<h at / = 100Hz can 
be measured for many binaries which do not have Si ? 2 
parallel to L. Though nominally an intrinsic parameter 
tied to the phase of the spins’ relative precession cycle, 
the parameter A<£> is limited by the range of precession 
dynamics allowed by the system at / re f. Second, the pos¬ 
terior in (J, £) is highly nongaussian, as shown in Figure 
[7] Finally, precession-induced modulation is known to 
break all degeneracies and enable strong constraints on 
many parameters EH E3 HU H3HS5]. As shown by 
the bottom panel of Figure [7] our ability to constrain 
the posterior distributions in (J, £) and in mass ratio is 
correlated. 

Figure [8] shows the posterior distributions for sources 
with different orientation Ojn relative to the line of 
sight. Similar to the change in morphology constraints 
with Ojn , the extent of the posterior distribution in the 
J£-plane increases (roughly) as sin Ojn- Nearly face-on 
sources have little precession-induced modulation in their 
GW signal, and therefore measurements cannot as reli¬ 
ably infer J, £. 

Finally, Figure [9] shows the posterior distributions for 
sources with different network amplitudes (SNR). Simi¬ 
lar to the change in morphology constraints with SNR, 
the extent of the posterior distribution in the J£-plane 
increases as the source SNR decreases. 

Comparing our results to the prior J, £ distribution 
shown in Figure [3] suggests the wide prior conven¬ 
tion has a significant impact on the posterior distribu¬ 
tion, “pulling” the posterior away from the region most 
strongly supported by the data. We anticipate that as- 
trophysically motivated priors which favor larger spins 
and more comparable mass ratios should produce tighter 
constraints on J, £ and morphology. 


C. Natural timescales as coordinates for 
double-spin measurements 

Rather than using the conserved constants (J, £) to 
characterize spin precession, we can use the two nat¬ 
ural precession frequencies ((f^/(27r)) , 1/r) introduced 
in Emm and reviewed in Appendix [B] Unlike the in¬ 
stantaneous precession vectors which appear in dSi ? 2 /dt, 
these expressions are independent of precession phase. 


Moreover, even for high-symmetry configurations like 
nearly aligned spins and post-Newtonian resonances, 
these timescales - corresponding to the rate of preces- 
sional modulation - are except for one special case well- 
defined and finite 0 though the associated amplitudes can 
be zero. For example, the frequency (Q z ) = a/r is 
nonzero for binaries with L,Si,S 2 all parallel, as can 
most easily be demonstrated in the special case of a sin¬ 
gle precessing spin PI- Similarly, even though no mod¬ 
ulation on the timescale takes place in post-Newtonian 
resonances, where spins and L remain coplanar, the spin 
precession timescale r is well-defined and nonzero [23|. 

As described in Appendix [B] due to coordinate am¬ 
biguity when J and L are parallel, the angle a can 
change by 2tt on two distinct surfaces in the J£-plane, 
but is otherwise continuous. The function r is continu¬ 
ous. Because the coordinate transformation from (J, £) 
to ((0, z ) /(27r), 1 /t) has two simple discontinuities, the 
connected posterior distributions in (J, £) seen above can 
map to two regions in ((fU) /(27r), 1/r) when the poste¬ 
rior distribution is sufficiently poorly constrained or close 
to these critical surfaces. More broadly, the transforma¬ 
tion between J,£ and Cl z ,r may not always be globally 
one-to-one: the transformation might lose some informa¬ 
tion. That said, because these timescales occur in the dy¬ 
namics and hence waveform, their imprint should appear 
in the posterior distribution, as with single-spin binaries 

m- 

Figure [lO] shows the posterior in these “observable” 
parameters: the precession frequencies (Q z ) and 1/r, 
evaluated at a GW frequency of 100Hz. First and fore¬ 
most, except for outliers associated with nearly aligned 
sources and templates, this figure suggests that, as a first 
approximation, GW measurements isolate a candidate 
source to a relatively narrow region of the ((fU),l/r) 
plane. These tight limits are only possible with suffi¬ 
cient leverage to identify both timescales in the modu¬ 
lated signal: when the viewing angle is small (Ojn — 0) 
or precession-induced modulation is small (L • J ~ 1, as 
with A<h = 7 r; see Figure [ 2 ]), then observations cannot 
pin down both timescales. Second, because (fl z ) = a/r 
is not a continuous function of (J, £), due to ambiguity 
in defining the precession angle a for some degenerate 
geometries, a posterior that is connected in the J£-plane 
can be disconnected in the ((Ct z ) , r) plane. That said, 
particularly for the narrowly confined A<h = 0 sources, 
the posterior in is connected. Third, for spin-orbit 
resonances, no relative spin precession occurs - all spins 
remain coplanar - so no modulation on the timescale r 
exists in the source spin dynamics; hence for all Ojn , no 
modulations on the timescale r can occur in the source 


6 As described in Gerosa et al. EH, for certain seperations a binary 
with both spins parallel to L, with the more massive “up” and 
the less massive “down”, is unstable to spin precession. When 
the instability occurs, a single point in the interior of the J, £ 
region has r — 00 . 
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FIG. 8. Constraining (J,£): dependence on Ojn • Posteriors from different simulations in the J,£ plane for fixed SNR 
and different Ojn , increasing from top (0jtv = 0) to bottom ( Ojn = 7t/2). Points are colored according to the £ values of the 
simulations; in this figure, all sources have SNR=20. 


waveform. We hypothesize that this lack of modulation, 
unique to our choice of sources, contributes to the rela¬ 
tive width of the posterior distributions in 1/r and (Cl z ), 
though we cannot verify our hypothesis with this sam¬ 
ple. Extrapolating from our sample, we further antici¬ 
pate that sources with precession timescales r compara¬ 
ble to or longer than the observationally accessible sig¬ 
nal duration T wave also cannot easily be distinguished 
from one another. Based on the limited sample avail¬ 
able in our investigation, coordinates that are constant 
on the precession timescale like (J, £) and ((fU), 1/r) are 
particularly valuable tools to extract robust statements 
about precessing binaries, isolating naturally correlated 
and tightly constrained parameters. Our methods gener¬ 
alize prior work using seperation-of-timescales to identify 
natural parameters and timescales for precessing black 
hole-neutron star binaries, and to relate those timescales 
to GW parameter estimation H3 eh eh]. We antici¬ 
pate these representations will be helpful when mining 
GW data from precessing BBHs for robust astrophysical 
statements. 


D. Biases, the maximum likelihood estimate, and 
results with noise 


Our study used doubly special sources. On the one 
hand, resonant binaries lie in a relatively small corner 
of the compact binary parameter space [22j 23 . On the 
other hand, the detector data we used was not a generic 
noise realization: for simplicity, we used exactly zero 
noise. As a result, as with posterior distributions of black 
hole spin and mass ratio for comparable-mass and highly 
spinning black holes, we generally cannot expect and do 
not observe the posterior confidence intervals in (J, £) or 
(Q Zl r) to be centered on the known source parameters: 
see e.g. Figs. [7] and [lOj Though our posterior distribu¬ 
tions had little support for the input binary’s true param¬ 
eters, our results are nonetheless well-converged and con¬ 
sistent with a single-best-fitting set of parameters close 
to the injected value. Figures [7] and [lO] both show a cross 
(x) at the location of the single sample with the max¬ 
imum likelihood. Unlike the mean posterior value, our 
maximum likelihood estimates are generally much closer 
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FIG. 9. Constraining (J,£): dependence on SNR. Posteriors from different simulations in the J^-plane for fixed 
Ojn = 7 t/4 and different SNR, increasing from top (SNR — 8) to bottom (SNR = 16). Points are colored according to the 
simulation’s £ value. 


to the injected value. Given how few posterior samples 
are available per run, we anticipate that observed dif¬ 
ferences between the maximum-likelihood estimates and 
the injected value are consistent with the limits of our 
finite sample size. 

Owing to our exceptional inputs and unusual poste¬ 
riors, to further demonstrate the stability of our results 
we repeated some of our analysis with a random noise 
realization. Figure [TT] shows our results. The posterior 
distributions remain highly nongaussian, with structures 
similar to those observed in the zero-noise study. 

E. Morphological classification: Data versus prior 

Unlike conclusions derived about the chirp masses and 
spins j4TM48] , for properties that may be difficult to mea¬ 
sure (such as the source morphology) the prior plays a 
surprisingly significant role. To quantify how much in¬ 
formation we learn from the data, we compute the in¬ 
formation gain (in bits) as the Kullback-Leibler (KL) 
divergence [66] between the morphological classification 


posterior and prior distributions: 

D K M\p*„m = X>w>>i°fe ( jSd ■ 

(7) 

where 6 = {A4>, £, Ojn, SNR}. The KL divergence has 
seen increasing application in physics mmol- If the pos¬ 
terior resembles the prior {p ~ p pr i or ), we have learned 
nothing from the data and the KL divergence is nearly 
0. To provide a sense of scale, for continuous gaussian 
distributions with the same mean but different standard 
deviations a ^ cr*, the KL divergence between these dis¬ 
tributions is Dkl — [ln(<r/<r*)] 2 / In 2; a KL divergence 
of order unity therefore implies a difference in mean by 
one standard deviation, or a difference in variance by a 
factor of order 2. Figure [12] shows this information gain 
as a function of the source parameters A4 >,£,#jtv. 

A comparison with Figure [4] reveals that, despite low 
posterior probabilities for A4> = i r sources, even in the 
best case {Ojn = tt/2 ) the peak amount of informa¬ 
tion we learn from the data is similar for the two source 
morphologies. The posterior probability depends signif- 
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FIG. 10. Constraining ((Q*), r) at 100 Hz. Posterior distributions for the A<f> = 0 sources (top left) and A<f> = tt sources 
(top right) in the ((Q z ) /(27t),1/t) plane, colored according to simulation £ value. Values of r and (Q z ) are evaluated using 
the binary configuration at the reference frequency / re f = 100 Hz. Solid (dashed) lines show 95% (67%) confidence intervals for 
each set of posterior samples from distinct simulations. Points are colored according to the £ values of the simulations which are 
marked by stars. Black crosses show, for each set of posterior samples, the maximum likelihood estimate. With few exceptions, 
the posterior associated with each AT = 0 injection is concentrated in a small range of (Q z ) and 1/r near the actual value. 
The results shown here are for Ojn — 7t/2, but similar results hold for 6jn > 7t/4. Bottom : Detail of the (Q z ) /(27t), 1/t plane, 
showing posteriors for simulations with £ > —0.4. 
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FIG. 11. Effect of noise on the posterior. As Figure [5] this figure shows the marginal posterior distribution in J, £ for 
sources with 6jn = 7t/4 and SNR=20, but including a random noise realization. Stars mark the actual values of J and £ from 
the simulations. Points and stars are and are colored according to the simulation’s £ value. 


icantly on the prior. Specifically, even though the data 
often strongly favors the relative rarity of [Eq. 
([6|] is responsible for the relatively poor classifications 
for A<£> = 7r sources, compared to A<£> = 0 sources. If we 
had adopted an astrophysically motivated prior, for ex¬ 
ample favoring comparable-mass binaries, we would have 
found a far more favorable result for the ability of GW 
measurements to distinguish between sources in distinct 
morphologies. We will address the impact of alternative 
priors in a subsequent study. 


F. Systematic uncertainty and conservation of £ 


The quantity £ is known to be conserved at 2PN or¬ 
der by the spin precession equations when the QM term 
m is included. This quantity may not be conserved at 
higher PN order, and it is known not to be conserved if 
the QM term is omitted. In this paper, we used the de¬ 
fault model for black hole spin precession implemented in 
lalsimulation, which explicitly omitted the QM term. 
This omission was not intentional - indeed, all previous 
analyses also omitted this term [4TH481 - but it provides 
an opportunity to assess the systematic error introduced 
when £ is not exactly conserved on the spin precession 
timescale, as could occur at higher PN order. 

Using the precession equations adopted in this analysis 

Haim, 


d(£LM 2 ) 

dt 


3L-Si x S 2 v 7 f 

-L • 

2r] 


[Si(i + q) — s 2(1 + 1 / o)] • 

(8) 


For configurations at or near a post-Newtonian reso¬ 
nance, £ will remain nearly constant, because L, Si and 
S 2 remain nearly coplanar. 

For this reason, even though the PN spin precession 
equations do not enforce it, for all source binaries £ is 
nearly constant throughout the evolution, rarely vary¬ 
ing by as much as 0.05. Compared to the typical sta¬ 
tistical measurement errors shown in Figures [7| [8] and [9j 
this systematic uncertainty is small for the signal ampli¬ 
tudes used in this study. Using a statistically significant 
subsample from each posterior distribution, we have also 
manually verified that £ is nearly conserved, occasionally 
oscillating but in all cases varying by much less than the 
statistical uncertainty in the posterior. Finally, we have 
repeated some of the simulations shown in Figure [7] with 
the QM term included and find quantitatively similar re¬ 
sults; see Figure [13] At least for the resonant sources 
used in this work, posterior parameter distributions are 
not sensitive to the inclusion of the QM term. 


IV. CONCLUSIONS 


In this work, we use the lalinference parameter es¬ 
timation code to infer the parameters of precessing BBHs 
from their GW signal. Contrary to some claims that the 
subdominant spin has little observationally accessible ef¬ 
fect m , we demonstrate by concrete example that the 
distinctive signature of uniquely two-spin physics can be 
inferred for the special but conservative set of sources we 
chose, featuring exactly resonant sources with mass ratio 
q = 0.8, Mtot = 13.5M 0 and maximal spins Specifically, 
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FIG. 12. Information gain. Kullback-Leibler (KL) divergence between the posterior and the prior distribution, measured 
in bits [Eq. 0 ]- Top: KL divergence as a function of £ and Ojn for fixed SNR = 20. Bottom: KL divergence as a function 
of £ and SNR for fixed Ojn — 7t/4. Structure in these plots reflects how the misalignment angles between the various angular 
momenta change versus £ and how this affects information gain: cf. Figures [2] and [4] As described in Appendix [C] the SNR=20 
results (red curve) were produced using a lower starting frequency than results for SNR< 20; this discrepancy is responsible 
for the slight difference in trend between the bottom left panel’s result for SNR=20 versus results for SNR< 20. 
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FIG. 13. Results using QM-corrected spin precession: After reproducing selected runs, we find similar posteriors in J, 
£ (top panel; compare to Figure [7]) and Q z ,r (bottom panel; compare to Figure 10 which adopts the same color scheme). 
Consistent with the corresponding figures, the top and bottom panel adopt different color schemes. 
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our analysis demonstrates that GW measurements can 
infer a subtle feature of the relative orientation of the 
two black hole spins: the binary’s precessional morphol¬ 
ogy. Each posterior parameter distribution can be sub¬ 
divided into three classes, identifying the relative proba¬ 
bility that the progenitor binary had one of three char¬ 
acteristic behaviors (librating about A4> = 0, librating 
about A4> = 7 r, or circulating). Using as sources the most 
extreme examples of the two librating cases - resonant 
binaries, which exhibit precession on only one of the two 
natural timescales - we confirmed the conclusion of [49] 
that resonant binaries can be correctly identified with the 
appropriate morphology. In other words, even though li¬ 
brating spins occupy a small part of parameter space, 
we show that these measurements reproducibly and cor¬ 
rectly identify the neighborhood associated to each reso¬ 
nance, except for finely tuned highly symmetric configu¬ 
rations. 

Because of computational cost our study focused on 
a single set of masses and spins. Further investigations 
with an astrophysically motivated distribution of masses, 
spins, and source distances are needed to assess how often 
and how reliably morphology can be classified in prac¬ 
tice. Our promising (if preliminary) results support the 
hypothesis that morphological classification could be an 
astrophysically interesting and robust diagnostic, as dis¬ 
tinct librating morphologies are both measurable and the 
natural consequence of distinct features of compact bi¬ 
nary evolution [[24] . 

Motivated by theoretical studies of precessing binaries 
[22] 23] , we present our results in a new coordinate sys¬ 
tem ((U 2 /(27t)) ,1/t), which better reflects the physical 
observables in the GW signal - precession-induced modu¬ 
lations in amplitude and phase - and respects precession 
physics, notably separation of timescales. Based on our 
limited study of special systems, this coordinate system 
seems to better characterize the dynamical information 
that GWs provide. 

Recently, Ref. ma used a small sample to investigate 
how well precessing BBH spins could be measured by Ad¬ 
vanced LIGO-scale instruments, using a small number 
of mass combinations, spin orientations, and distances 
(3x2x3), but several viewing angles Ojn- Like previ¬ 
ous studies, the result of Vitale et al. uni showed that 
posterior distributions depend sensitively on viewing an¬ 
gle, being best constrained when seen edge-on. Addi¬ 
tionally, they investigated how well subdominant spins 
could be measured for precessing binaries, claiming on 
the basis of a limited sample that the relative angle A<f> 
between the projection of the two spins into the orbital 
plane cannot be measured. Our results provide a concrete 
counterexample. That said, Ref. US used very differ¬ 
ent initial configurations with no resonant sources; they 
adopted suboptimal diagnostics for the second spin’s ori¬ 
entation (i.e., the time-dependent azimuthal angle A4>, 
like Fig. [5| instead of morphology as in Fig. [4] or J, £ 
point as in Fig. [ 7 ]); and they had too few independent 
initial conditions to robustly explore the full parameter 


space. Specifically, their sample included a handful of 
sources: three mass pairs, only one of which involved 
unequal-mass BH-BH binaries (1OM 0 + 5M 0 ); fixed BH 
spin magnitudes; and only two relative spin orientations, 
only one of which was not in He (i.e., not circulating). 
For their single noncirculating source - a (1OM 0 ,5M 0 ) 
pair of rapidly (xi = X 2 = 0-9) spinning black holes with 
initial spins Si oriented by 45° and S 2 oriented by 135° 
from the line of sight and coplanar with L - these authors 
do note that for Ojn ~ 7t/2 the posterior distribution of 
A4> is significantly different from the prior. As noted 
above, further study with an astrophysically motivated 
distribution of masses, spins, source distances, and ori¬ 
entations is required to assess how well GWs can identify 
the source morphology of generic nonresonant sources. 

Our analysis employed inspiral-only waveforms which 
lack the coalescence and ringdown signals present in real 
BBH merger signals. Their unphysical termination con¬ 
ditions are known to introduce convent ion-dependent ar¬ 
tifacts into parameter estimation, with increasing impact 
as the total binary mass increases HD [72]. A detailed 
discussion of waveform termination conditions is beyond 
the scope of this paper. That said, based on the rela¬ 
tively high termination frequencies shown in Table [TT| we 
anticipate that waveform termination conditions do not 
dominate the differences we observe. 


ACKNOWLEDGMENTS 

D.T. is partially supported by the National Sci¬ 
ence Foundation through awards PHY-1067985, PHY- 
1404139, PHY-1055103 and PHY-1307020. D.T. is grate¬ 
ful for the support and hospitality of V. Kalogera’s 
group and the Center for Interdisciplinary Exploration 
and Research in Astrophysics (CIERA) at Northwest¬ 
ern University, where this project was conceived. R.O’S. 
is supported by NSF Award PHY-1505629. D.G. is 
supported by the UK STFC and the Isaac Newton 
Studentship of the University of Cambridge. E.B. is 
supported by NSF CAREER Grant PHY-1055103 and 
by FCT contract IF/00797/2014/CP1214/CT0012 under 
the IF2014 Programme. M.K. is supported by Alfred 
P. Sloan Foundation grant FG-2015-65299. T.B.L. ac¬ 
knowledges NSF award PHY-1307020. U.S. is supported 
by FP7-PEOPLE-2011-CIG Grant No. 293412, FP7- 
PEOPLE-2011-IRSES Grant No.295189, H2020-MSCA- 
RISE-2015 Grant No. StronGrHEP-690904, H2020 
ERC Consolidator Grant Agreement No. MaGRaTh- 
646597, SDSC and TACC through XSEDE Grant 
No. PHY-090003 by the NSF, Finis Terrae through 
Grant No. ICTS-CESGA-249, STFC Roller Grant No. 
ST/L000636/1 and DiRAC’s Cosmos Shared Memory 
system through BIS Grant No. ST/J005673/1 and STFC 
Grant Nos. ST/H008586/1, ST/K00333X/1. Computa¬ 
tional resources were provided by the Northwestern Uni¬ 
versity Grail cluster (CIERA) through NSF MRI award 
PHY-1126812, by the Atlas cluster at AEI Hannover, 


20 


supported by the Max Planck Institute and by the Nemo 
at cluster through NSF-092340. 


Appendix A: Coordinate conventions for precessing 

spins 


For reference and clarity, in this appendix we summa¬ 
rize the different coordinate conventions used to describe 
precessing spins: the radiation frame, in which angular 
momenta are parametrized relative to the line of sight; 
the “system frame” [39] 08], in which angular momenta 
are parametrized relative to the total angular momen¬ 
tum J; the spin parameters { 6 lS\ ? @ls 2 ? frequently 
used in the resonant-locking literature, in which the spin 
angular momenta are parametrized relative to the or¬ 
bital angular momentum L; and the (J, £, S) and (J, £, p') 
coordinate sytems introduced in previous work [22, [23] . 
Unless otherwise noted, all parameters are specified at 
/ref = 100 Hz. 

The radiation frame m expresses the orbital and spin 
angular momenta in polar coordinates relative to the di¬ 
rection of GW propagation — N, with the zero of az¬ 
imuthal angle set by the plane spanned by L, — N. 

In the system frame [39] , the angular momenta of the 
binary are characterized in a hierarchical fashion. The 
angles Ojn and 0j are polar coordinates of the total an¬ 
gular momentum J direction relative to the direction of 
propagation —N; one angle 0 jl characterizes the rela¬ 
tive orientation of L on its precession cone around the 
total angular momentum, equivalently constraining the 
direction of the total spin S = Si + S2; and three an¬ 
gles (t \, £ 2 , 0 i 2 ) characterize the relative orientation of 
the two component spins relative to the orbital angular 
momentum. For example, t\ = cos -1 Si • L and, if x and 
y are unit vectors that form a right-handed orthonormal 
frame with L, e 2 ^ 12 = A x [(x + iy) • S 2 ]/[(x + iy) • Si] 
where A is a real amplitude chosen so the norm of the 
right-hand side is unity. 

The resonant-locking literature uses parameters 
#i, 02 > A<I> to characterize the relative orientation of the 
two spins relative to L. These parameters are identical 
to the tilt and relative phase angles (t\ = 0 \ = OlSh 
t<i = 62 = 0ls 2 i an d A<f> = 0i2) used in the system frame 
currently favored for lalinference [39]. 

Finally, recent work |[22J, 23, provides an explicit solu¬ 
tion for the relative two-spin dynamics, expressing spin 
evolution in the three-dimensional relative spin space 
(e.g., 0 ls 1 , 0 ls 2 i A<F) in yet another triplet of coordinates 
J, £, 5 on the three-dimensional space of all relative spin 
orientations at fixed L, Si, S 2 , and q. Gerosa et al. [23] 
also introduce an alternative coordinate system «/,£,<//. 
In these two systems, the first coordinate, the magni¬ 
tude J of the total angular momentum, is a conserved 
constant on precession timescales; the projected effective 
spin £ [cf. Eq. §] is conserved on all timescales by the 


orbit-averaged 2PN spin precession equation^] the co¬ 
ordinate S is the magnitude of the total spin; and the 
coordinate p' is an angle characterizing the degree to 
which the plane spanned by Si, S 2 is rotated relative to 
the plane spanned by S, L. Specifically, the angle p' is de¬ 
fined by decomposing the spins relative to S = Si + S 2 ; 
a vector perpendicular to S, L; and the remaining vector 
of an orthonormal triad [Fig. 1 in [23]]: 

2 < ai > 

In terms of these coordinates, we can redundantly express 
the spin using the spin magnitude S and the relative 
angle p' between the plane spanned by the two spins and 
the plane spanned by S and L: 


51 = AS + 6(cos p'X + sin p'Y) (A2a) 

5 2 = (1 - A) S - b(cosp'X + sin p'Y) (A2b) 


s^ + sl-sl 

2S 2 

b 2 = S\ ~ A 2 S 2 


1 -a= s 2±A^1 (A2c) 


2S 2 


(A2d) 


Either the phase angle p' or the spin S can be eliminated 
from these expressions; see Gerosa et al. [2T for details. 


Appendix B: Spin precession formulae 


For each value of mi,m2 ,Xi,X2 5 ^,/0 Refs. [22l [23] 
define two quantities that characterize precession: a and 
r. The timescale r is the duration of one precession cycle 
of the spins relative to L, defined by 


r = 



dS 

\dS/dt\ 


(Bl) 


where S is the magnitude of the total spin (cf. Appendix 
\Al and where ( dS/dt) 2 is a polynomial in S with simple 
roots at the turning points S = S±. Because l/\dS/dt\ 
has an integrable singularity at both limits S±, r is finite 
and continuous everywhere, except under certain condi¬ 
tions at the single point in the (J, £) plane associated with 
the larger spin aligned and the smaller spin antialigned 
with L [61 . By contrast, the angle a is the azimuthal 
angle accumulated during the precession of L around J. 
The azimuthal angle 0 jl of L around J is is therefore not 
well defined at times when L and J are parallel. As an in¬ 
tegral of d0jL/dt, a inherits this ambiguity. Away from 
this surface, the function a is well-defined, continuous, 
and can be evaluated using the definite integral 


a = 





dS 


\dS/dt\ 


(B2) 


7 Using the complete 2PN spin precession equations, without orbit 
averaging, £ varies on the orbital timescale. 








21 


where the rational function Q z is given by Eq. (19) in 
[22] . This discontinuity in a arises from an integrable 
singularity in Cl z ; it implies a discontinuity in a/r, and 
it is responsible for the E-shapes seen in Figure [TO] 
Because the J£-plane and (Q z ) = a/r play a significant 
role in our presentation, we provide an explicit formula 
to determine where these discontinuities occur. The con¬ 
ditions on (J, £) needed for colinearity follow from Eqs. 
(6) and (14) in |23], after substituting suitable values of 
S. For clarity, we also derive them from straightforward 
vector algebra. Because colinearity of L and J implies all 
angular momenta are coplanar, without loss of generality 
we use the polar angles 0js 1 0js 2 £= [0,7r] of Si ? 2 relative 
to J to characterize the two spins. In these coordinates, 
colinearity of L and J requires 


J = sL + Si cos 6j Sl + S 2 cos 0js 2 (B3) 

0 = Si sin 0 JSl + S 2 sin 0js 2 (B4) 

§ 

^ = M t miXl COS ° JSl + 1712X2 COS ° JS ^ ( B5 ) 


where s = dzl corresponds to two possible choices (align¬ 
ment or antialignment) between the orbital and total an¬ 
gular momentum: J = 5 L 

Keeping in mind that 0js 1 ,0js 2 G [0, i r], we find with¬ 
out loss of generality that these conditions imply one 
quadratic equation and one system of two linear equa¬ 
tions on the variables cosO js 1 , cos 0js 2 ‘ 

Si( 1 - cos 2 6^) = Sf( 1 - cos 9jg 2 ) (B6) 


' J - sL 


r 9 91 

mf 777-2 

Xi cos 0js 1 



777-1 /M m 2 /M 

_X2 COS djs 2 _ 


Solving the linear equations for cos 0js 1 , cos 0js 2 and 
then substituting into the quadratic constraint generally 
leads to a quadratic form in (J, £). For this problem, 
however, the coefficient of £ 2 cancels, leaving a linear ex¬ 
pression for £ as a function of J. The solution [compare 
e.g. to Eq. (13) in [23]] reads 

, = V - L *) 2 

* S 2r/M 2 (J — Ls) 

(1 + q) 2 ( J ~ Ls) 2 - (1 - q 2 )(Sf - S%) , 

= s - W^lW 2 -' ( 8) 


Using this solution, we find that the angles 0js 1 ^0js 2 
vary with J according to 


cos 0js 1 = 


COS 0JS 2 = 


[(J-£ a) 2 + ( ^-S 2 2)] 

2Si( J — Ls) 

[(J - Ls) 2 - (Sf - £%)] 
2 S 2 (J~Ls) 


(B9) 

(BIO) 


Not all values of J correspond to realizable angles for 
0js 1 ,0js 2 - For geometric reasons - and as can be im¬ 
mediately verified by direct substitution of appropriate 
choices for J, namely J = L =b S\ =b S 2 - the curve with 
5 = 1 passes through all four aligned-spin configurations 



FIG. 14. Alignment between L and J and disconti¬ 
nuities in a in the J£-plane. As in Figure [3] the thick 
black lines show the allowed region in the J£-plane for our 
fiducial parameters, and the dashed lines show J c o,ctt(£)- The 
orange curves show the locus of points J, £ where L and J are 
parallel and thus where a is discontinuous, as determined by 
solutions to Eq. (B8) that correspond to physically realizable 
spin orientations. 


IcosfljsJ = |cos0js 2 | = 1. Conversely, for s = 1, a 
region containing J = L is always excluded. In the lan¬ 
guage of [23], the root 5 = 1 corresponds to = \J—L\ 
and the root 5 = — 1 corresponds to S ma ^ «= J + L. For 
example, as illustrated by Figure [14] for our fiducial pa¬ 
rameters, one of these curves (orange in the figure) con¬ 
nects the point cos 6 js x = cos 0 js 2 = 1 to the “kink” in 
the J c o(£) curve corresponding to 0js 1 = 7 t,0js 2 = 0, 
the “up/down” instability point where r = 00 [61]. The 
other curve joins the remaining two spin-aligned con¬ 
figurations and is nearly equal to the minimum value 
J m in(0 allowed at each £ over the corresponding range 
of £. This latter discontinuity is responsible for the sig¬ 
nificant change in Q z at £ = £t seen in Table [IT] and for 
the V-shape feature seen in Figure [To] 


Appendix C: Technical details and caveats 

Small technical details about waveform generation and 
sampling can have a dramatic impact on final results. In 
this appendix, we provide more exhaustive details about 
the specific calculations we perform. 

All signals are sampled with an interval At = 1/2048s, 
corresponding to a Nyquist frequency f^ yq = 1024Hz. 
Our signal and templates start at 2/ or b = 10Hz. [For 
sources with SNR < 20, we used 2/ or b = 20Hz.] Each bi- 
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A<f> 

£ 

J / M 2 

/end (Hz) 

A0jl/27T 

Ojl 

0ls x 

0ls 2 

(Q z )/2tt[Hz] 

l/r [Hz] 

0 

- 1.0 

0.39 

281 

3.50 

0.00 

7T 

7T 

- 

- 

0 

- 0.9 

0.49 

303 

3.98 

0.46 

2.63 

2.79 

2.13 

2.16 

0 

- 0.8 

0.58 

321 

4.47 

0.55 

2.41 

2.63 

2.46 

2.07 

0 

- 0.7 

0.66 

345 

4.83 

0.58 

2.23 

2.51 

2.73 

1.99 

0 

- 0.6 

0.72 

375 

5.20 

0.59 

2.09 

2.39 

2.95 

1.92 

0 

- 0.5 

0.78 

403 

5.59 

0.59 

1.95 

2.29 

3.14 

1.84 

0 

- 0.4 

0.84 

438 

5.94 

0.58 

1.83 

2.19 

3.30 

1.77 

0 

- 0.3 

0.89 

483 

6.25 

0.56 

1.71 

2.09 

3.44 

1.71 

0 

- 0.2 

0.94 

533 

6.58 

0.55 

1.60 

1.99 

3.55 

1.65 

0 

- 0.1 

0.99 

608 

6.88 

0.52 

1.50 

1.90 

3.65 

1.59 

0 

0.0 

1.03 

687 

7.18 

0.50 

1.39 

1.80 

3.74 

1.54 

0 

0.1 

1.08 

763 

7.47 

0.47 

1.29 

1.69 

3.81 

1.49 

0 

0.2 

1.12 

812 

7.73 

0.45 

1.19 

1.59 

3.87 

1.44 

0 

0.3 

1.16 

812 

7.99 

0.42 

1.09 

1.48 

3.91 

1.40 

0 

0.4 

1.20 

834 

8.22 

0.39 

0.99 

1.36 

3.95 

1.36 

0 

0.5 

1.23 

823 

8.48 

0.35 

0.88 

1.23 

3.97 

1.33 

0 

0.6 

1.27 

825 

8.67 

0.31 

0.77 

1.10 

3.99 

1.30 

0 

0.7 

1.30 

831 

8.90 

0.27 

0.66 

0.95 

3.99 

1.26 

0 

0.8 

1.34 

865 

9.11 

0.22 

0.53 

0.77 

3.99 

1.23 

0 

0.9 

1.37 

857 

9.29 

0.16 

0.37 

0.54 

3.97 

1.21 

0 

1.0 

1.40 

843 

9.47 

0.00 

0.00 

0.00 

- 

- 

7T 

- 1.0 

0.39 

281 

6.29 

0.00 

7T 

7T 

- 

- 

7r 

- 0.9 

0.44 

299 

6.44 

0.04 

2.84 

2.55 

- 

- 

7T 

- 0.8 

0.48 

317 

6.62 

0.06 

2.74 

2.27 

3.92 

1.92 

7T 

- 0.7 

0.53 

336 

6.83 

0.08 

2.69 

2.04 

3.91 

1.80 

7T 

- 0.6 

0.58 

359 

7.05 

0.09 

2.67 

1.81 

3.92 

1.71 

7T 

- 0.5 

0.62 

384 

7.30 

0.10 

2.68 

1.58 

3.94 

1.65 

7 T 

- 0.4 

0.67 

413 

7.57 

0.10 

2.72 

1.33 

3.96 

1.61 

7T 

- 0.3 

0.71 

444 

7.92 

0.09 

2.78 

1.05 

3.99 

1.60 

7T 

- 0.2 

0.75 

481 

8.29 

0.07 

2.89 

0.71 

4.02 

1.60 

7T 

- 0.1 

0.79 

518 

4.51 

0.06 

2.92 

0.10 

2.47 

1.59 

7T 

- 0.0 

0.86 

569 

4.80 

0.16 

2.43 

0.32 

2.60 

1.48 

7T 

0.1 

0.92 

636 

5.12 

0.19 

2.14 

0.44 

2.71 

1.38 

7T 

0.2 

0.98 

712 

5.44 

0.20 

1.91 

0.53 

2.79 

1.28 

7T 

0.3 

1.04 

755 

5.72 

0.19 

1.69 

0.60 

2.85 

1.20 

7T 

0.4 

1.10 

799 

5.96 

0.17 

1.49 

0.65 

2.89 

1.13 

7T 

0.5 

1.15 

819 

6.23 

0.15 

1.30 

0.67 

2.91 

1.08 

7T 

0.6 

1.20 

834 

6.46 

0.13 

1.11 

0.66 

2.90 

1.06 

7T 

0.7 

1.25 

850 

6.65 

0.10 

0.92 

0.62 

2.89 

1.07 

7T 

0.8 

1.30 

864 

6.80 

0.08 

0.72 

0.53 

2.86 

1.09 

7T 

0.9 

1.35 

879 

6.91 

0.05 

0.49 

0.39 

2.82 

1.13 

7T 

1.0 

1.40 

843 

6.98 

0.00 

0.00 

0.00 

- 

- 


TABLE II . Derived source properties. For each dynamically distinct binary used as a source of GWs, this table reports 
on properties of the dynamics, including: /end, the termination frequency of the GW signal; A (/)jl, the number of precession 
cycles of L around J between the starting and ending frequencies; Ojl , Ols x , Ols 2 , characterizing the relative orientations of the 
angular momenta J,L,Si,S 2 ; J , the magnitude of the total angular momentum; and (Q z ) /(27r) and 1/r, the characteristic 
precession frequencies. The last five columns are evaluated at the reference frequency, 100Hz. Discontinuities in {Q z } occur as 
described in Appendix [5] 


nary is evolved until terminated at the minimum stable 
circular orbit or when the orbital frequency begins to de¬ 
crease with time, whichever comes first. For the masses 
and spins studied here, this termination frequency is sig¬ 
nificantly smaller than the Nyquist frequency, but signif¬ 
icantly higher than the frequency at which most power 
accumulates; see, e.g., Fig 4 in [49] , 

Our signals were analyzed in a 128s data window, far 
longer than the signal. [For sources with SNR < 20, we 
used 16s.] As described elsewhere, rather than use the 
full 15-dimensional likelihood, we explicitly marginalized 


it over time and orbital phase at each step in the MCMC 
[48] . Using selected examples, we have confirmed that 
our results are unchanged if this marginalization is not 
used. As seen in Table [TT] particularly for £ ~ —1, some 
of our signals were short, being terminated at compara¬ 
tively small frequencies due to the breakdown of the post- 
Newtonian expansion we employ to generate them. That 
said, not only is relatively little signal power associated 
with / > 300Hz, but also no structure in the posterior 
correlates tightly with the termination frequencies listed 
in Table [TT} A detailed study of the impact of termination 
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conditions on precessing parameter estimation is beyond 
the scope of our work. 

With the exception of Figure mi in the text we use 
an exactly zero-noise realization. The likelihood of a 
given set of GW detector data depends on the detec¬ 
tor noise in two ways B3 EE]: through the specific de¬ 
tector data being analyzed, and through the probabil¬ 
ity of any given noise realization. Each instrument’s 
noise n(t) is assumed to be stationary and gaussian 
Markov process, characterized by a noise power spec¬ 
trum ( n(u:')*n(uj )) = Sh(u)S(uj — o/)/2; using the noise 


power spectrum, we can evaluate the probability of any 
noise realization n(t). To synthesize a unique data set 
to be analyzed for each set of intrinsic parameters A, 
we assume that the data d(t) in each instrument - gen¬ 
erally containing both signal and noise, or of the form 
d(t) = h(t\X) + n(t) for each instrument — is given ex¬ 
actly by d{t) = h{t) = F+h+(t\\) + F x h x (t\X), where 
F +?x are the detector response functions and h+ :X (t\X) 
are the two linear polarizations of of the GW. Just like 
evaluating a normal distribution at the mean, this arbi¬ 
trary choice eliminates ambiguity in subsequent Bayesian 
calculations of the posterior distribution. 
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